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Abstract 

We summarize our recent result for a splitting function for small x evolution which 
includes resummed small x logarithms deduced from the leading order BFKL equation with 
the inclusion of running coupling effects. We compare this improved splitting function with 
alternative approaches. 
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In recent years the theory of scaling violations for deep inelastic structure functions at small 
x has attracted considerable interest, prompted by the experimental information coming 
from HERA. New effects beyond the low-order perturbative approximation to anomalous 
dimensions or splitting functions should become important at small- a;. However, no major 
deviation of the data from a standard next-to-leading order perturbative treatment of 
scaling violations has been found. By now the origin of this situation has been mostly 
understood [1,2]. 

The BFKL kernel x(ct s ,M) has been computed to next-to-leading accuracy (NLO): 

X (M,a s ) = a sX o(M) + a 2 sX i(M) + .... (1) 

The problem is how to use the information contained in xo an d Xi in order to improve 
the splitting function derived from the perturbative leading singlet anomalous dimension 
function 7(a s , N) which is known up to NLO in a s : 

7 (iV,a s ) =a s7o (iV) + a 2 sll (N) + (2) 

in such a way that the improved splitting function remains a good approximation down to 
small values of x. This can be accomplished [1] by exploiting the fact that the solutions 
of the BFKL and GLAP equations coincide at leading twist if their respective evolution 
kernels are related by a "duality" relation. In the fixed coupling limit the duality relation 
is simply given by: 

X (l(N,a s ),a s ) = N. (3) 

The splitting function then will contain all relative corrections of order (a s log l/x) n , de- 
rived from xo(M), and of order a s (a s log l/x) n : derived from xi(M). 

The early wisdom on how to implement the information from xo was completely 
shaken by the computation of xi, which showed that the naive expansion for the improved 
anomalous dimension had a hopelessly bad behaviour. However, as a consequence of 
the physical requirement of momentum conservation, this problem is cured if the small- 
x resummation is combined with the standard resummation of collinear singularities, by 
constructing a 'double-leading' perturbative expansion. 

However, higher order corrections to the kernel qualitatively change the asymp- 
totic small-x behaviour of structure functions by changing it from x~ x ° to 
x~ x = x~ x °e Ax t w x~ x °[l + AA£ + ....], with: 

A = A + AA, A = a s Xo(|) = " s c , AA = a 2 s xi(\) H . (4) 

The computed correction Xi{\) to the leading result is quite large, and this suggests 
that the correct asymptotic exponent is not reliably determined by the two known terms. 
Therefore in our previous work we have treated A as a parameter to be fitted from the 
data. Good agreement with the HERA data was found, but at the price of having to 
sharply fine-tune this parameter. 

In our recent work [1] we have shown that also this problem can be solved by a 
full treatment of running coupling corrections. Indeed, it has been known for some time 
that running coupling effects can be included perturbatively order by order at small x by 
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adding effective subleading Axi contributions to the BFKL kernels XiiXii — However, 
these additional terms turn out to have singularities at M = 1/2, which correspond to an 
enhancement by powers of Inl/x of the associated splitting functions, which may offset 
the perturbative suppression by powers of a s . In our approach, we have shown that these 
contributions can be resummed to all orders at the level of splitting functions, in a way 
compatible with factorization and a smooth behaviour in the small- £ limit. 

Based on these results, we now know the way the information contained in xo(M) 
should be used in order to construct a better first approximation for the improved anoma- 
lous dimension. Indeed, we find that, once running coupling effects are properly included 
in the improved anomalous dimension, the asymptotic behavior near x = is much soft- 
ened with respect to the naive Lipatov exponent. Hence, the corresponding dramatic rise 
of structure functions at small x, which is phenomenologically ruled out, is replaced by a 
milder rise. The effective exponent corresponding to this rise turns out to be in surpris- 
ingly good agreement with the fine-tuned value of A indicated by the fit, which in turn led 
to a splitting function which closely followed the NLO GLAP in the region of the data. 
This suggests that additional higher order corrections are small, and that a leading-order 
approximation based on the standard BFKL kernel xo is phenomenologically viable. 

We discuss now explicitly our proposed improved splitting function. Assuming that 
one only knows 70 (AT), 71 (AT) and xo(M), the improved anomalous dimension has the 
following expression: 

7 f L (a s ,iV) = [a sl0 (N) + a 2 sll (N) + ls (jf) - 2^] + 

1 as (5) 

+ 7a(c ,« s , N) - 2 + y^r s i N ~ a sCo] + i/3oa s (l + -j^c ) - mom. sub. 

The first line on the right-hand side, within square brackets, is the double-leading expres- 
sion for the improved anomalous dimension at this level of accuracy, made up of the NLO 
perturbative term a s ^o(N) + a^ji(N) plus the power series of terms (a s /N) n in 7s (^f ), 
obtained from xo using eq. (3), with subtraction of the order a s term to avoid double 
counting (ca = n c = 3). In the second line, the "Airy" anomalous dimension 7a (co, « s , N) 
contains the running coupling resummation, and the remaining terms subtract the contri- 
butions to 7a(co, « s , N) which are already included in 7 S , 70 and 71. The Airy anomalous 
dimension 7,4(00, a s , N) is the exact solution of the running coupling BFKL equation cor- 
responding to a quadratic appoximation of xo near M = |: Xo ~ [co + \kq(M — \) 2 }. 
Finally "mom. sub." is a subleading subtraction that ensures momentum conservation 
1I {a s ,N= 1) = 0. 

We now summarize the properties of the improved anomalous dimension in this ap- 
proximation. In the limit a s — > with arbitrary N, 7/(a s ,A r ) reduces to a s ^o(N) + 
Qig7i(A r ) + O(ctg). For a s — > with a s /N fixed, 7/(a s ,A r ) reduces to 7dl-lo = 
«s7o(A r ) + 7 S (jf) — 2 t^, i-e. the leading term of the double-leading expansion. Thus 
the Airy term is subleading in both limits. In spite of this, its role is very significant be- 
cause of the singularity structure of the different terms in eq. (5). In fact, 70 (N) has a pole 
at N = 0, 7 S has a branch cut at = a s co, and 7^4 has a pole at = Ao < a s co, where 
A^o is the position of the rightmost zero of the Airy function. The importance of the Airy 
term is that the square root term subtracted from 7^4 cancels, within the relevant accuracy, 
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Figure 1: The improved splitting function corresponding to 7/(a s , N) eq. (5) with a s = 
0.2 (dot-dashed), compared with those from the DL-LO approximation (dashed) and 
GLAP NLO (solid). 



the branch cut of 7 S at N = a s co and replaces the corresponding asymptotic behaviour 
at small x with the much softer one from 7^4. Note that the quadratic approximation is 
sufficient to give the correct asymptotic behaviour up to terms which are of higher order 
in comparison to those included in the double-leading expression in eq. (5). 

While we refer the interested reader to ref. [1] for further details, we show here two new 
figures for the splitting function which illustrate the success of this approach and compare 
it both to our older results and to those of ref. [2]. In figure 1 we show (for a s = 0.2) the 
singlet splitting function obtained from eq. (5) compared with the NLO GLAP kernel and 
with the DL-LO approximation, which displays the sharp small- a; rise characteristic of the 
BFKL resummation. In the region of the HERA data, our improved splitting function, 
with no free parameters, closely follows the NLO GLAP evolution with a behaviour at small 
x which is much softer than that of BFKL. It is interesting to note that the agreement 
between GLAP and resummed results is significantly improved by the inclusion in eq. (5) 
of 71 and the corresponding double-counting subtraction. This improvement was already 
shown in ref. [1] in the anomalous dimension, and is even more apparent in the splitting 
function shown here. 

In figure 2 the improved splitting function is compared with those of our earlier ap- 
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Figure 2: The improved splitting function corresponding to -fi(a s , N) eq. (5) with a s = 
0.2 (dot-dashed), compared with GLAP NLO (solid), with the result of our old approach, 
where the value of A = 0.21 was fine tuned to best fit the HERA data (dashed), and 
with the splitting function of ref. [2] (dotted). Note that the vertical scale is rather 
larger than that in figure 1. 



proach to this problem and with that recently obtained in ref. [2]. It is seen that the 
correct inclusion of running coupling corrections is by itself sufficient to produce the soft- 
ening of the behaviour at small x required in order to reproduce the data. Thus it is no 
longer necessary to introduce A by hand as a free parameter. In fact, the fine-tuned value 
A ~ 0.21 (for a s ~ 0.2) obtained from fitting the data is in remarkable agreement with the 
value iVo = 0.2112 determined from the Airy anomalous dimension with a s = 0.2. This 
improved splitting function is so good that subleading corrections to it are presumably 
very small. 

In the same figure, we also show the splitting function from ref. [2]. In comparison 
to the approach adopted there, we share the general physical framework, but there are a 
number of differences. Specifically, the authors of ref. [2] determine the asymptotic small- 
er behaviour by making some assumptions on the way to implement a symmetrisation 
M — > (1 — M) of the BFKL kernel x( CK s, M), which we prefer not to do because of the 
ambiguities it introduces. Also, our resummed curve in figure 2 includes the effects of 71 but 
not fully those of xi- This is because we find that after running coupling resummation, 
the impact of the left-over terms in xi is comparable to the resummation ambiguities. 
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By contrast, their curve includes the effects of xi but not those of 71. Furthermore, in 
ref. [2] the running coupling evolution equation is solved explicitly but numerically to the 
relevant order. Instead, through the Airy expansion, we obtain an analytic expression 
which holds to the same accuracy. This opens up the possibility to fit the data and match 
to usual evolution equations in a straightforward way. The resummed splitting functions 
shown in figure 2 all rise asymptotically with approximately the same power, but differ in 
the preasymptotic region. These differences can be taken as a measure of the remaining 
resummation uncertainties. 
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